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PEST/WALD/PINIT: ADAPTIVE PSYCHOPHYSICS TESTING PACKAGE 


DECUS Program Library Write-up 


DECUS NO. 8-520 


/THIS IS A PACKAGE CONSISTING OF THREE PROGRAM LISTINGS - 
/WALD, PEST, AND PINIT. THESE PROGRAMS IMPLEMENT THE 
/PROCEDURES DESCRIBED IN "PEST: EFFEICINT ESTIMATES ON 
/PROBABILITY FUNCTIONS" (TAYLOR AND CREELMAN, JOURNAL OF THE 
/ACOUSTICAL SOCIETY OF AMERICA, 1967, VOLUME 41, PAGES 
/782-787). THE PROGRAMS ARE WRITTEN FOR A PDP-8/S COMPUTER 
/BUT WILL WORK ON ANY 8-SERIES MACHINE. 

/THE USE OF THE PROGRAMS IS AS FOLLOWS: ON EACH OF A NUMBER 
/OF EXPERIMENTAL TRIALS, THE COMPUTER PRESENTS A STIMULUS TO 
/A SUBJECT AND ASKS HIM TO MAKE SOME KIND OF DECISION. THIS 
/DECISION CAN BE OBJECTIVELY SCORED AS CORRECT OR INCORRECT. 
/AND THE FUNCTION OF THE PEST PACKAGE IS TO DETERMINF THAT 
/STIMULUS LEVEL WHICH PRODUCES A DESIRED TARGET PROPORTION 
/OF CORRECT RESPONSES. 


/THIS DOCUMENTATION DEALS WITH PROGRAMMING, RATHER THAN WITH 
/THEORETICAL, ASPECTS OF THE USE OF PEST. THE PACKAGE IS AN 
/ENTIRELY SELF-CONTAINED SET OF SUBROUTINES FOR PROCESSING 
/TRIAL-BY-TRIAL RESULTS OF AN EXPERIMENT IN CONJUNCTION WITH 
/A 4-WORD WALD PARAMETER LIST AND A 7-WORD PEST HISTORY LIST 
/FOR EACH EXPERIMENT BEING RUN. THE WALD PARAMETER LIST 
/SPECIFIES THE TARGET PERCENTAGE AND THE DEVIATION FROM THAT 
/PERCENTAGE AT WHICH A DECISION TO CHANGE LEVELS IS MADE, 
/ALONG WITH THE CURRENT STATUS OF THE STATISTIC USED TO 
/DECIDE IF THE CURRENT PROPORTION OF TRIALS CORRECT IS 
/WITHIN THE BOUNDS ESTABLISHED IN THE PARAMETER LIST. 

/THE PEST HISTORY LIST CONTAINS THE CURRENT LEVEL OF THE 
/STIMULUS PARAMETER BEING ADJUSTED AND INFORMATION ABOUT 
/THE RANGE OF ALLOWABLE ADJUSTMENTS, IN ADDITION TO SOME 
/INFORMATION ABOUT THE MOST RECENTLY-MADE ADJUSTMENTS. 

/FROM THE STANDPOINT OF THE PEST PACKAGE, AN EXPERIMENT IS 
/WHATEVER REFERENCES A GIVEN PAIR OF WALD AND PEST PARAMETER 
/LISTS. THUS, ANY NUMBER OF CONCURRENT EXPERIMENTS CAN BE 
/RUN WITH A SINGLE 1-PAGE PEST PACKAGE, PLUS 11 CORE 
/LOCATIONS PER EXPERIMENT. 


/THE WALD PROGRAM IS ENTERED WITH A 1 OR 0 FOR A SUCCESS OR 
/FAILURE, RESPECTIVELY. IT RETURNS A 0 IF NO DECISION HAS 
/YET BEEN MADE, A +1 IF IT HAS BEEN DECIDED THAT THE 
/PROPORTION OF SUCCESSES IS TOO HIGH, AND A -1 IF THE 
/PROPORTION IS TOO LOW. THESE NON-0 OUTPUTS CAN THEN BE USED 
/AS INPUTS TO PEST ITSELF, WHICH DECIDES WHAT NEW STIMULUS 
/VALUE TO USE AND INSTALLS THAT VALUE INTO THE PEST HISTORY 
/LIST AND ALSO RETURNS IT IN THE ACCUMULATOR. IF THE CHANGE 
/IN LEVEL RESULTED IN ACCUMULATOR OVERFLOW OR UNDERFLOW, THE 
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/LINK IS RETURNED SET, OTHERWISE CLEAR. THE ONLY OTHER 
/EXTERNALLY-CALLED ROUTINE IN THE PEST PACKAGE IS THE PEST 
/INITIALIZER, PINIT, WHICH TAKES A LIST OF FOUR CONSTANTS 
/AND COPIES IT IN ANOTHER LOCATION WITH THREE 0S APPENDED, 
/SO THAT PEST CAN OPERATE ON A LIST OF VARIABLES WITHOUT 
/DESTROYING THE LIST OF INITIALIZING CONSTANTS. 

/PEST IS NORMALLY RUN IN "MOUSE" MODE (MININUM OVERSHOOT AND 
/UNDERSHOOT SEQUENTIAL ESTIMATION). IN THIS MODE, WHEN 
/A DECISION WOULD BE MADE TO MAKE A STEP SMALLER THAN THE 
/SPECIFIED MINIMUM, THE PROGRAM MAKES A NON-STANDARD 
/SUBROUTINE RETURN TO INDICATE THAT A FINAL LEVEL HAS BEEN 
/REACHED. A SIMPLE PATCH ALLOWS THE PROGRAM TO BE RUN IN 
/"RAT" MODE (RAPID ADAPTIVE TRACKING), IN WHICH A STEP 
/SMALLER THAN THE ALLOWED MINIMUM IS SIMPLY REPLACED WITH 
/A MINIMUM STEP AND THE PEST PROGRAM HAS NO STOPPING RULE. 
/IN THIS MODE, HOWEVER, THE LINK CANNOT 3E USED FOR OVERFLOW 
/INFORMATION. DETAILS OF. THE PROGRAM PATCH FOLLOW PEST IN 
/THE PROGRAM LISTING. 

PAUSE 
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*1600 / /NOMINAL — THIS IS FULLY 

/ / /RELOCATABLE TO ANY CORE PAGE 


&LUM TOUE&f 08 *' ™ * SUCCESS ' ACC:0 * MI1 -URE. 
/ JMS WALD 

/ WALD0 / /ADDRESS OF WALD PARAMETER LIST 

/ (PROGRAM RETURNS HERE, USUALLY TO A TEST FOR A 

/ NON-ZERO ACCUMULATOR, IN WHICH CASE PEST IS CALLED) 

/THE PROGRAM RETURN HAS ACC=+1 FOR TOO MANY SUCCESSES -7 
/FOR TOO FEW, AND 0 FOR WITHIN BOUNDS. IF THE DEVIATION 
/EQUALS THE LIMIT EXACTLY, IT IS CONSIDERED OUT OF BOUNDS. 
/THE WALD PARAMETER LIST CONSISTS OF FOUR ENTRIES, THE FIRST 
/THREE OF WHICH ARE CONSTANTS USED BUT NOT CHANGED BY WALD. 
/THE LAST THE ONLY LOCATION NEEDED AS A WALD HISTORY SPAC* 7 . 
/THE FOUR LOCATIONS ARE, IN ORDER: 

/WALD0, NUMERATOR OF EXPECTED FRACTION CORRECT 
/ COMMON DENOMINATOR OF EXPECTATION AND LIMIT 

/ NUMERATOR OF DEVIATION LIMIT 

/ SPACE FOR WALD HISTORY 

/THE WALD HISTORY LOCATION SHOULD BE INITIALIZED BY SETTING 
/IT TO 0 FOR WALD TO CONSIDER THE FIRST TRIAL IT RECEIVES. 
/OR TO -(DEVIATION LIMIT NUMERATOR), TO IGNORE THAT TRIAL 
/ALL THREE NUMBERS MUST BE POSITIVE INTEGERS IN THE RANGE 
/FROM 1 TO 3777 OCTAL (1 TO 2047 DECIMAL), AND FURTHER. 

/THE SUM OF THE LAST TWO (DEVIATION LIMIT NUMERATOR + COMMON 
/DENOMINATOR) MUST NOT EXCEED 3777 OCTAL (2047 DECIMAL). 

/FOR EXAMPLE, A WALD PARAMETER LIST TO AIM FOR A TARGET 
/PERCENTAGE OF 75% WITH A DEVIATION LIMIT OF 1.5 
/TRIALS COULD HAVE ENTRIES: 

/WALD0, 3 /NUMERATOR OF EXPECTATION 

/ 4 - /COMMON DENOMINATOR 

/ 6 /NUMERATOR OF DEVIATION LIMIT 

/ 0 /LOCATION FOR WALD HISTORY 

/FRACTIONS NEED NOT BE EXPRESSED IN LOWEST TERMS 


/ONCE WALD HAS BEEN RUNNING, IT WILL RESET THE HISTORY 
/LOCATION TO -(DEVIATION LIMIT NUMERATOR) WHENEVER THE 
/DEVIATION LIMIT HAS BEEN EXCEEDED, AND WILL IGNORE THE 
/NEXT TRIAL. 


1600 

0000 WALD, 

0 

1S01 

3315 

DCA WPRES 

1602 

1600 

TAD I WALD 

1603 

4343 

JMS WPSET7 
/ / 

1604 

2200 

ISZ WALD 

1605 


/DECIDE IF THIS 

1766 

TAD I WHISTP 

1606 

1765 

TAD I WDNUMP 

1607 

7640 

SZA CLA / 


/ / 

/ / 


/STORE RESULTS OF TRIAL 
/GET ADDRESS OF PARAMETER LIST 
/SETUP A LIST OF SEVEN CONSECUTIVE 
/POINTERS TO SEVEN LOCATIONS 
/CORRECT THE RETURN ADDRESS 

TRIAL SHOULD BE IGNORED 

/ADD HISTORY TO DEVIATION NUMERATOR 
/IF HISTORY = - DEV LIMIT, THAT IS 
/THE INDICATOR THAT THIS WAS 
/THE FIRST TRIAL OF A SET, AND 
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/ 

/ 

/IS THEREFORE TO BE IGNORED. 

1S10 

5213 

JMP 

WPDOIT 

/IF NOT EQUAL, GO PROCESS THE TR1 

1611 

3766 

DCA 

I WHISTP 

/IF EQUAL, RESET HISTORY TO 0 

1612 

5600 

JMP 

I WALD 

/EXIT WITH A 0 FOR NO CHANGE 



/UPDATE THE ACCUMULATED DEVIATION FROM EXPECTATIC 

1613 

1315 

WPDOIT, TAD 

WPRES 

• /GET THE RESULT AGAIN 

1S14 

7640 

SZA 

CLA / 

/AND, IF IT INDICATES A SUCCESS 

1615 

1764 

TAD 

I WDENP 

/GET THE DENOMINATOR 

1616 

7041 

CIA 

/ 

./NEGATE THAT OR 0 

1617 

1763 

TAD 

I WENUMP 

/ALWAYS ADD THE EXPECTATION 



/ 

/ 

/NUMERATOR 

1620 

1766 

TAD 

I WHISTP 

/ADD TO CURRENT HISTORY 

1621 

3766 

DCA 

I WHISTP 

/AND STORE THAT HISTORY 



/DECIDE IF THE 

DEVIATION EXCEEDS THE LIMIT 

1622 

1766 

TAD 

I WHISTP 

/NOW INSPECT IT AGAIN 

1623 

7500 

SMA 



1624 

704 1 

CIA 

/ 

/MAKE IT NEGATIVE 

1625 

1765 

TAD 

I WDNUMP 

/ADD TO DEVIATION LIMIT 

1625 

7740 

SMA 

SZA CLA 

/RESULT STAY NEGATIVE OR 0? 

1627 

5600 

JMP 

I WALD 

/IF NOT, WITHIN LIMITS, SO EXIT W] 



/ 

/ 

/ACCUMULATOR SET TO 0 



/IF 

A DECISION 

HAS BEEN REACHED, RESET THE HISTORV 



/LOCATION TO CAUSE THE NEXT TRIAL TO BE IGNORED. 



/ANE 

> EXIT WITH 

+1 OR -1, AS APPROPRIATE 

1630 

1765 

TAD 

I WDNUMP 

/FIRST, 

1631 

7041 

CIA 



1632 

3766 

DCA 

I WHISTP 

/RESET THE HISTORY LOCATION 

1633 

7040 

CMA 

/ 

/NOW SETUP A -1 

1634 

1315 

TAD 

WPRES 

/AND, IF TRIAL WAS A SUCCESS, 

1635 

1315 

TAD 

WPRES 

/CHANGE TO +1 

1636 

5600 

JMP 

I WALD 

/AND EXIT 


PAUSE 



1637 0000 

1640 3315 

1641 1637 

1642 4343 


1643 2237 


1644 

1645 

1646 

1647 


1650 

1651 

1652 

1653 

1654 


1767 

1315 

7640 

5256 


1770 

3771 

3770 
1764 
71 10 


/ENTER PEST WITH A +1 IN THE ACCUMULATOR FOR TOO MANY 
/SUCCESSES, A -1 FOR TOO MANY FAILURES (THE SAME AS THE 
/OUTPUT OF WALD). DO NOT EVEN CALL PEST IF NO CHANGE IS TO 

FOLLOWS: ™ E STIHUU,S - ™ E CALLING SMUMCE 

/ JMS PEST 

/ P .£o T oT«t/ /ADDRESS OF PEST PARAMETER LIST 

/ JMP STOP /PROGRAM RETURNS HERE IF STOPPING 

f f t /RULE CRITERION IS MET 

/ JMS SETLEV /PROGRAM RETURNS HERE WITH NEW LEVEL 

' ' / /IN THE ACCUMULATOR IF TRIALS ARE TO 

/ / / /CONTINUE • 


/THE PEST PARAMETER LIST CONSISTS OF THESE SEVEN LOCATIONS* 
/PEST0, <1 CURRENT STIMULUS LEVEL, NO RESTRICTIONS 

(2) CU RR ENT STEP SIZE, POSITIVE IN RANGE 1-3777 OCTAL 

(3) MINIMUM STEP SIZE, POSITIVE IN RANGE 1-3777 OCTAL 

(4) MAXIMUM STEP SIZE, POSITIVE IN RANGE 1-3777 OCTAL 

(5) -l OR +1, AS LAST STEP RESULTED FROM TOO MANY 

/ FAILURES OR SUCCESSES, INITIALIZED TO 0 

(6) -l OR 0, AS LAST STEP WAS DOUBLING OR NOT. 

/ INITIALIZED TO 0 

(7) COMSECUTIVE STEPS IN SAME DIRECTION, SUBTRACT 1 

/ IF REVERSAL FOLLOWED NON-DOUBLED STEP, 

/ SUBTRACT 2 IF IT FOLLOWED A DOUBLING 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


PEST, 0 


DCA WPRES 
TAD I PEST 
JMS WPSET7 
/ / 

/ . / 

ISZ PEST 
/ / 


/STORE ACCUMULATOR CONTENTS 
/GET ADDRESS OF PEST PARAMETER LIST 
/SETUP LIST OF SEVEN CONSECUTIVE 
/POINTERS TO SEVEN CONSECUTIVE 
/LOCATIONS 

/SETUP PEST TO POINT TO RETURN IF 
/STOPPING CRITERION MET 


/DECIDE IF THIS IS A REVERSAL OF DIRECTION 
TAD I PDIRPT 

/ADD NEW DIRECTION TO LAST DIRECTION 
/DO THEY CANCEL? 

/IF NOT, CONTINUE IN SAME DIRECTION 


TAD WPRES 
SZA CLA / 
JMP PSAME 


PDIFF, 


/PROCESS REVERSAL OF DIRECTION 


1655 5267 


TAD I PDUBPT 
DCA I PSMCTP 
/ / 

DCA I PDUBPT 
PDFLAT, TAD I PSSPT 
CLL RAR / 

/ 

/ 

/ 

/ 


/ 

/ 

/ 

/ 


JMP PSTNSS 


/GET DOUBLING POINTER, 0 OR -1 
/PLACE IN COUNTER OF CONSECUTIVE 
/STEPS IN SAME DIRECTION 
/STORE 0 FOR NOT DOUBLED STEP 
/GET CURRENT STEP SIZE 
/DIVIDE BY TWO AND TRUNCATE 
/NOTE: IF A LOW-ORDER BIT IS LOST BY 
/THIS TRUNCATION, IT CANNOT BE 
/RESTORED BY A SUBSEQUENT STEP-SIZE 
/DOUBLING 

/GO STORE THE SIZE 
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/PROCESS A CONTINUATION IN THE SAME DIRECTION 


1656 

2771 

PSA ME, 

ISZ 

I PSMCTP 


/INCREASE THE COUNT OF CONSECUTIVE 




/ 

/ 


/STEPS IN THIS DIRECTION 




/ 

/ 


/WHEN THE COUNT REACHES 2, THAT 




/ 

/ 


/MEANS AT LEAST 3 STEPS IN THIS 




/ 

/ 


/DIRECTION, SINCE THE FIRST STEP 




/ 

/ 


/ONLY INITIALIZES BUT DOES NOT 




/ 

/ 


/INCREMENT THE COUNT. IF THE SKIP 




/ 

/ 


/OCCURS, THE RESULT IS 0, WHICH IS 




/ 

/ 


/NOT STRICTLY POSITIVE, AND SO IT 




/ 

/ 


/DOES NOT MATTER THAT THE -1 ABOUT 




/ 

/ 


/TO BE CREATED IS MISSED. 

1657 

7240 


CLA 

CMA / 


/CREATE A -1 

1660 

1771 


TAD 

I PSMCTP 


/ADD THE COUNT 

1661 

7750 


SPA 

SNA CLA 


/STRICTLY POSITIVE YET? 

1662 

5266 


JMP 

PSTNSS-1 


/IF NOT, DON’T DOUBLE 

1663 

7040 


CMA 




1664 

3770 


DCA 

I PDUBPT 


/IF SO, STORE -1 FOR DOUBLING 

1665 

1764 


TAD 

I PSSPT 


/FETCH THE STEP SIZE v 

1666 

1764 


TAD 

I PSSPT 


/EITHER ONCE OR TWICE ^ 

1667 

3764 

PSTNSS, 

DCA 

I PSSPT 


/STORE NEW STEP SIZE 




/IF 

THE STEP 

IS 

TOO LARGE, CUT IT IN HALF 

1670 

1766 


TAD 

I PMAXPT 


/FETCH MAXIMUM SIZE 

1671 

7140 


CLL 

CMA / 


/CLEAR LINK AND CREATE -(MAXIMUM+I) 




/ 

/ 


/NOW STEP SIZE EQUAL TO MAXIMUM IS 




/ 

/ 


/OK; STEP SIZE ''GREATER TRIPS LINK 

1672 

1764 


TAD 

I PSSPT 


/ADD CURRENT STEP SIZE 

1673 

7630 


SZL 

CLA / 


/A NONZERO LINK MEANS ONE OF TWO 




/ 

/ 


/THINGS, EITHER THE STEP SIZE IS 




/ 

/ 


/LARGER THAN 3777, OR LARGER THAN 




/ • 

/ 


/THE MAXIMUM ALLOWABLE, OR BOTH 

1674 

5253 


JMP 

PDFLAT 


/SO CUT IT DOWN AGAIN 




/IF 

THE STEP 

IS 

\ 

TOO SMALL, STOP CHANGING LEVELS 

1675 

1765 


TAD 

I PMINPT 


/NOW TRY THE MINIMUM 

1676 

7041 


CIA 




1677 

1764 


TAD 

I PSSPT 


• 

1700 

7710 

PRATMD, 

SPA 

CLA / 


/POSITIVE OR EQUAL MEANS STEP SIZE ^ 




/ 

/ 


/IS NOT TOO SMALL 

1701 

5637 


JMP 

I PEST 


/TOO SMALL MEANS RETURN TO LOCATION 




/ 

/ 


/MEANING STOPPING RULE IS SATISFIED 




/THE STEP IS 

ALLOWABLE, SO TAKE THAT STEP AND RETURN 

1702 

2237 


ISZ 

PEST 


/NOW MAKE RETURN ADDRESS POINT TO 




/ 

/ 


/NORMAL LEVEL-CHANGE LOCATION 

1703 

7100 


CLL 

/ 


/PREPARE TO RECEIVE OVERFLOW 

1704 

1315 


TAD 

WPRES 



1705 

3767 


DCA 

I PDIRPT 


/STORE LATEST DIRECTION 

1706 

1764 


TAD 

I PSSPT 


/GET STEP SIZE 

1707 

2315 


ISZ 

WPRES 


/TEST DIRECTION 

1710 

7041 


CIA 

/ 


/NEGATE FOR TOO MANY SUCCESSES 

1711 

1763 


TAD 

I PLEVPT 


/ADD TO CURRENT LEVEL 

1712 

3763 


DCA 

I PLEVPT 


/STORE 
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1713 

1714 


1763 

5637 


TAD I PLEVPT 
JMP I PEST 


/GET BACK AGAIN 
/EXIT WITH IT 


/TO RUN THIS PROGRAM IN 
/PATCHES: 

/PRATMD,SMA CLA / 

/ / / 

/ JMP .+3 / 

/ TAD I PMINPT 

/ / / 

/ DCA I PSSPT 

/ 

/AND USE A SHORTER SUBROUTINE CALL: 

/ JMS PEST 

/ PEST0 / /ADDRESS OF PARAMETER LIST 

/ (SUBROUTINE RETURNS HERE WITH NEW LEVEL) 

PAUSE 


"RAT" MODE, MAKE THE FOLLOWING 

/POSITIVE OR EQUAL MEANS STEP SIZE 

/IS NOT TOO SMALL 

/SO GO CHANGE THE LEVEL 

/IF STEP WOULD BE TOO SMALL, REPLACE 

/IT WITH THE MINIMUM 
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/THE ROUTINE PINIT IS CALLED AS FOLLOWS: 




/ 

JMS 

PINIT 




/ 

PPERM / 

/ADDRESS OF STARTING VALUES OF FIRST 



/ 

/ 

/ 

/FOUR ITEMS OF PEST LIST 



/ 

PEST0 / 

/ADDRESS OF A BLOCK OF 7 LOCATIONS 



/ 

/ 

/ 

/INTO WHICH THE FOUR ITEMS ABOVE 



/ 

/ 

/ 

/CAN BE COPIED, FOLLOWED BY 3 ZEROS 

1715 

0000 

PINIT, 

0 



1716 

1715 


TAD 

I PINIT 

/GET ADDRESS OF PERMANENT LIST 

1717 

3237 


DCA 

PINITP 

/STORE IN A POINTER 

1720 

2315 


ISZ 

PINIT 


1721 

1715 


TAD 

I PINIT 

/GET ADDRESS OF WORKING LIST 

1722 

4343 


JMS 

WPSET7 

/SETUP 7 POINTERS 

1723 

2315 


ISZ 

PINIT 

/CREATE THE RETURN ADDRESS 

1724 

1637 


TAD 

I PINITP 

/GET FIRST PERMANENT ITEM 

1725 

3763 


DCA 

I PLEVPT 

/STORE IN WORKING LIST 

1726 

2237 


ISZ 

PINITP 


1727 

1637 


TAD 

I PINITP 


1730 

3764 


DCA 

I PSSPT 

(ft 

1731 

2237 


ISZ 

PINITP 

W 

1732 

1637 


TAD 

I PINITP 


1733 

3765 


DCA 

I PMINPT 


1734 

2237 


ISZ 

PINITP 


1735 

1637 


TAD 

I PINITP 


1736 

3766 


DCA 

I PMAXPT 

/SAME FOR SECOND, THIRD, FOURTH 

1737 

3767 


DCA 

I PDIRPT 


1740 

3770 


DCA 

I PDUBPT 


1741 

3771 


DCA 

I PSMCTP 

/NOW 3 ZEROS 

1742 

5715 


JMP 

I PINIT 

/EXIT 



/ROUTINE TO 

SETUP SEVEN 

CONSECUTIVE POINTERS TO SEVEN 



/CONSECUTIVE CORE LOCATIONS 

1743 

0000 

WPSET7, 

0 



1744 

3363 


DCA 

WENUMP 

/STORE FIRST POINTER, GIVEN ON ENTRY 

1745 

1372 


TAD 

WPM6 


1746 

3373 


DCA 

WPCNTR 

/GET A -6 AND STORE AS COUNTER 

1747 

1374 


TAD 

WPTAD 

/GET A FETCHING INSTRUCTION ^ 

1750 

3353 


DCA 

WPTADP 

/PLACE IT IN THIS ROUTINE 

1751 

1375 

• 

TAD 

WPDCA 

/AND A DEPOSITING INSTRUCTION 

1752 

3355 


DCA 

WPDCAP 

/PLACE IT IN ITS SPOT 

1753 

7000 

WPTADP, 

NOP 

/ 

/FETCH A POINTER 

1754 

7001 


IAC 

/ 

/INCREMENT THE POINTER 

1755 

7000 

WPDCAP, 

NOP 

/ 

/DEPOSIT IT IN THE NEXT CELL AS A 

. 



/ 

/ 

/POINTER TO THE NEXT LOCATION 

1756 

2353 


ISZ 

WPTADP 


1757 

2355 


ISZ 

WPDCAP 

/INCREMENT BOTH INSTRUCTIONS 

1760 

2373 


ISZ 

WPCNTR 

/DO THIS LOOP SIX TIMES 

1761 

5353 


JMP 

.-6 / 

/DO IT AGAIN IF NECESSARY 

1762 

5743 


JMP 

I WPSET7 

/EXIT IF DONE 


/THE SEVEN BASIC POINTERS: 


WENUMPr. / /WALD EXPECTATION NUMERATOR 

1763 0030 PLEVPT, 0 / /PEST CURRENT LEVEL 


8 




* 

WDENPr. / 

/ 

1764 

0000 

PSSPT, 0 

/ 



WDNUMPr. 

/ 

1765 

0000 

PMINPT, 0 

/ 



WHISTPs. 

/ 

1766 

0000 

PMAXPT, 0 

/ 

1767 

0000 

PDIRPT, 0 

/ 

1770 

0000 

PDUBPT, 0 

/ 

1771 

0000 

PSMCTP, 0 

/ 



WPRES=PINIT 

/ 

1772 

7772 

WPM6, -6 


1773 

0000 

WPCNTR, 0 


1774 

1363 

WPTAD, TAD 

WENUMP 



/ 

/ 

1775 

3364 

WPDCA, DCA 

WENUMP+1 



/ 

/ 



/ 

/ 



PINITPrPEST 



/WALD COMMON DENOMINATOR 
/PEST CURRENT STEP SIZE 

/WALD DEVIATION LIMIT NUMERATOR 
/PEST MINIMUM STEP SIZE 

/WALD HISTORY 

/PEST MAXIMUM STEP SIZE 

/PEST DIRECTION 

/PEST DOUBLING INDICATOR 

/PEST CONSECTIVE STEPS INDICATOR 

/WHICH WILL BE FREE THEN 


/INSTRUCTION TO FETCH FIRST OF THE 
/SEVEN POINTERS, FOR INSERTION ABOVE 
/INSTRUCTION TO DEPOSIT INTO SECOND 
/OF THE SEVEN POINTERS, FOR 
/INSERTION ABOVE 




PD FLAT 

1653 

PDIFF 

1650 

PDIRPT 

1767 

PDU3PT 

1770 

PEST 

1637 

PI NIT 

1715 

PINITP 

1637 

PLEVPT 

1763 

PMAXPT 

1766 

PMINPT 

1765 

PRATMD 

1703 

PSA ME 

1656 

PSMCTP 

1771 

PSSPT 

1764 

PSTNSS 

1667 

WALD 

1600 

WDENP 

1764 

WDNUMP 

1765 

WENUMP 

1763 

WHISTP 

1766 

WPCNTR 

1773 

WPDCA 

1775 

WPDCAP 

1755 

WPDOIT 

1613 

WPM6 

1772 

WPRES 

1715 

WPSET7 

1743 

WPTAD 

1774 

WPTADP 

1753 



/DEMONSTRATOR PACKAGE FOR WALD AND PEST 
/ /HOWARD KAPLAN - 1971 

/THIS PACKAGE CONSISTS OF ROUTINES THAT CAN BE USED TO 
/TEST, DEMONSTRATE, AND SIMULATE THE ACTIONS OF WALD AND 
/PEST. THE WALD TEST ROUTINE GETS INPUT FROM ANY SUBROUTINE 
/THAT GENERATES 1 AND 0 (CORRECT AND INCORRECT) RESPONSES. 
/THIS PACKAGE INCLUDES TWO SUCH ROUTINES. ONE SIMPLY ACCEPTS 
/"l" OR ”0” INPUTS FROM THE KEYBOARD. THE OTHER GENERATES 
/A +1 WITH ANY SPECIFIED PROBABILITY, AFTER A USER-WRITTEN 
/SUBROUTINE HAS CONVERTED THE CURRENT STIMULUS LEVEL TO A 
/PROBABILITY, USING ANY DESIRED PSYCHOMETRIC FUNCTION. 

/WHEN THE WALD TEST REACHES ITS DECISION, IT PRINTS ITS 
/RESULT AND EITHER RECYCLES FOR A DECISION ON ANOTHER STRING 
/OF INPUTS OR ELSE PASSES CONTROL TO THE PEST PACKAGE FOR 
/A CHANGE OF LEVEL, AGAIN UNDER THE CONTROL OF ONE POINTER. 
/IT IS ALSO POSSIBLE TO SUPPLY THE PEST PACKAGE WITH +1 
/OR -1 INPUTS FROM ANY OTHER SUBROUTINE, SUCH AS THE SIMPLE 
/TELETYPE INPUT ROUTINE SUPPLIED. 

/THE DEMONSTRATOR PACKAGE USES A COPY OF THE PEST PACKAGE 
/LOCATED AT 1600, AND ITSELF OCCUPIES SPACE FROM 2000- 
/2777. WHILE THE PACKAGE AS WRITTEN DOES NOT USE THE 
/FLOATING POINT SYSTEM, THIS PLACEMENT INTO THE LOW END OF 
/CORE ALLOWS ANY USERS TO EXTEND ITS SIMULATION CAPABILITIES 
/TO USES REQUIRING FLOATING-POINT ARITHMETIC FOR SUCH 
/USES AS CALCULATING NEW SUCCESS PROBABILITIES. THE ONE-PAGE 
/PEST PACKAGE ITSELF MAY BE RELOCATED, WITHOUT CHANGE, TO 
/ANY CORE PAGE. 

/THIS PROGRAM IS DESIGNED TO BE USED WITH ODT LOW, FOR SUCH 
/USES AS CHANGING TO NEW STARTING PARAMETERS AND NEW 
/POINTERS BETWEEN ROUTINES. IN PARTICULAR, TO RUN THIS 
/PROGRAM, THE FOLLOWING MUST BE DONE: 

/l. PLACE THE APPROPRIATE CONSTANTS INTO THE WALD PARAMETER 
/LIST TO SPECIFY THE TARGET PERCENTAGE AND DEVIATION LIMIT. 
/THIS LIST BEGINS AT LOCATION WALD0, AND, AS AN EXAMPLE, 
/CONTAINS THE VALUES NEEDED FOR A TARGET PROPORTION OF .75 
/AND A DEVIATION LIMIT OF 1.0. 

/2. PLACE THE APPROPRITE CONSTANTS INTO THE PEST PARAMETER 
/LIST TO SPECIFY THE INITIAL (STIMULUS) LEVEL, INITIAL STEP 
/SIZE, MINIMUM STEP SIZE, AND MAXIMUM STEP SIZE. THIS LIST 
/BEGINS AT LOCATION PEST0, AND, AS AN EXAMPLE, CONTAINS 
/THE VALUES NEEDED FOR AN INITIAL LEVEL OF 100, AN INITIAL 
/STEP SIZE OF 16, A MAXIMUM STEP SIZE OF 32, AND A MINIMUM 
/STEP SIZE OF 2. 

/3. PLACE THE APPROPRIATE POINTERS INTO LOCATIONS WINPUT AND 
/PINPUT, TO SPECIFY WHERE THE INPUTS TO WALD AND TO PEST 
/WILL BE GENERATED. AS LOADED, WALD GETS ITS INPUT FROM THE 
/TELETYPE, AND PEST GETS ITS INPUT FROM THE RESULT OF THE 
/WALD TEST. 

/4. PLACE THE APPROPRIATE POINTERS INTO LOCATIONS WNEXT AND 
/PNEXT, TO SPECIFY WHETHER EACH ROUTINE SIMPLY LOOPS 
/BACK TO ITSELF OR ALTERNATES WITH THE OTHER WHEN IT IS 
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/DONE. ONCE THE PROGRAM IS STARTED, THE ROUTINES RE-ENTER AT 
/WLOOP AND PLOOP, RESPECTIVELY. THE PROGRAM AS A WHOLE 
/STARTS AT WALDGO OR AT PESTGO. 

/5. IF SUBROUTINE PROB1 IS BEING USED, THEN PLACE INTO 
/LOCATION CONVRP t. POINTER TO A SUBROUTINE WHICH WILL 
/CONVERT THE CURRENT LEVEL INTO A PROBABILITY FOR ROUTINE 
/PROB1 TO USE. FOR EXAMPLE, IF A CLASSICAL PSYCHOMETRIC 
/FUNCTION IS BEING USED (CUMULATIVE NORMAL), THEN THE 
/SUBROUTINE SHOULD CONVERT THE CURRENT LEVEL TO THE 
/AREA OF THE CURVE TO THE LEFT OF THE CURRENT LEVEL, IN AN 
/APPROPRIATE SCALE. THE NUMERATOR AND DENOMINATOR OF THIS 
/PROBABILITY SHOULD BE LEFT IN LOCATIONS PNUM AND PDEN. 

/A VERY ROUGH APPROXIMATION TO SUCH A FUNCTION IS PROVIDED 
/BY THE ROUTINE CALLED "CURVE", AT 2600. 

/FOR FURTHER INSTRUCTIONS, CONSULT THE PROGRAM LISTING. 

/A DEMONSTRATION OF THE USES OF THIS PROGRAM IS INCLUDED 
/AFTER THE LISTING. 

PAUSE 
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★ 2003 



/START 

HERE 

TO BEGIN 

2300 

4734 

WALDGO, 

JMS 

I FIXALP 

2001 

4753 

WLOOP, 

JMS 

I RESETP 




/ 

/ 




/ 

/ 

2002 

4753 

WLINE, 

JMS 

I RESETP 

2003 

1335 


TAD 

M71 

2334 

3336 


DCA 

LCOUNT 

2005 

2336 

WAGAIN, 

ISZ 

LCOUNT 

2006 

7413 


SKP 


2007 

5202 


JMP 

WLINE 

2013 

4775 


JMS 

I WINPUT 

2011 

3340 


DCA 

RESULT 

2012 

1340 


TAD 

RESULT 

2013 

1354 


TAD 

C263 

2014 

4752 


JMS 

I TYPEP 

2015 

1340 


TAD 

RESULT 

20 IS 

4747 


JMS 

I WALDP 


W 2017 2204 WALDO / 


2023 

7450 


SNA 

/ 

2021 

5205 


JMP 

WA GAI N 

2022 

3340 


DCA 

RESULT 

2023 

1340 


TAD 

RESULT 

2024 

704 1 


CIA 


2025 

1341 


TAD 

C254 

2026 

4752 


JMS 

I TYPEP 

2027 

5773 


JMP 

I WNEXT 



/START 

HERE 

TO BEGIN 

2033 

4734 

PESTGO, 

JMS 

I FIXALP 

2331 

3340 


DCA 

RESULT 

2332 

4753 

PPRINT, 

JMS 

I RESETP 

2033 

1340 


TAD 

RESULT 

2034 

7440 


SZA 

/ 

2035 

5245 


JMP 

PRINTC 

2036 

1342 


TAD 

M7 

2037 

3340 


DCA 

RESULT 

2040 

1337 


TAD 

SPACE 

2041 

4752 


JMS 

I TYPEP 

2042 

2343 


ISZ 

RESULT 

2043 

5240 


JMP 

.-3 / 

2044 

5251 


JMP 

PRINTL 

2045 

1341 

PRINTC, 

TAD 

C254 




/ 

/ 




/ 

/ 

2046 

4752 


• JMS 

I TYPEP 

2047 

1364 


TAD 

PSTEPW 

2050 

4746 


JMS 

I SPRINP 

2051 

1363 

PRINTL, 

TAD 

PLEVW 

2052 

4746 


JMS 

I SPRINP 

2053 

1370 


TAD 

PDUBW 

2054 

7650 


SNA 

CLA / 

2055 

1343 


TAD 

C56 / 


WITH A WALD TEST 

/DO INITIALIZATIONS 

/DO AN EXTRA CARRAIGE RETURN BEFORE 

/THE FIRST LINE OF A NEW WALD TEST 

/RECORD 

• /GO TO A NEW LINE 

/ALLOW 70 CHARACTERS PER LINE 
./LINE FULL YET? 

/IF SO, START A NEW ONE 
/GET A 0 OR 1 
/SAVE TEMPORARILY 


/ECHO RESULT 

/GO TO WALD ROUTINE 

/WITH THE RESULT 

/USING THIS PARAMETER LIST 

/NON-0 RESULT ON RETURN? 

/IF NOT, GET MORE INPUT 

/SAVE RESULT 


/PRINT "+" FOR +1, FOR -1 

/AND DO WHATEVER COMES NEXT 

PEST WITH PRINTOUT OF INITIAL LEVELS 
/DO INITIALIZATIONS 
/NO RESULT KNOWN YET 
/START A NEW LINE 
/GET RESULT 
/IS THERE ONE? 

/IF SO, PRINT CHANGE OF LEVEL 


/IF NOT, TYPE 7 SPACES 

/THEN GO PRINT STARTING LEVEL 

/PRINT "+" FOR RESULT OF -1, FOR 

/RESULT OF +1, TO GIVE PROPER 

/SIGN TO AMOUNT OF LEVEL CHANGE 

/TYPE SIGN 

/GET CHANGE 

/type in decimal, then 2 spaces 

/PRINT CURRENT LEVEL 
/IS DOUBLING POINTER 0? 

/IF SO, PRINT "N" FOR NOT DOUBLED 
/THIS CHANGES SPACE TO "N" 
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2056 

1337 


TAD 

SPACE 

/IF NOT, PRINT 3LANK 

2057 

4752 


jms 

I TYPEP 


2060 

1344 


TAD 

LETRD 

/PRINT ”D” IN EITHER CASE 

2061 

4752 


JMS 

I TYPEP 


2062 

1337 


TAD 

SPACE 


2063 

4752 


JMS 

I TYPEP 


2064 

1337 


TAD 

SPACE 


2065 

4752 


JMS 

I TYPEP 

/AND 2 SPACES 

2066 

1371 


TAD 

PCONW 

/NOW GET COUNT OF CONSECUTIVE 




/ 

/ 

/STEPS IN THIS DIRECTION 

2067 

7700 


SMA 

CLA 


2073 

5277 


JMP 

POSCT 

/GO AWAY IF POSITIVE 

2071 

1345 

NEGCT, 

TAD 

MINUS 


2072 

4752 


JMS 

I TYPEP 

/PRINT SIGN 

2073 

1371 


TAD 

PCONW 


2074 

7041 


CIA 



2075 

4746 


JMS 

I SPRINP 

/AND PRINT NEGATIVE OF COUNT 

2076 

5303 


JMP 

PFIX 

/AND THAT'S IT 

2077 

1337 

POSCT, 

TAD 

SPACE 


2100 

4752 


JMS 

I TYPEP 

/TYPE A SPACE 

2101 

1371 


TAD 

PCONW 


2102 

4746 


JMS 

I SPRINP 

/AND THE CONSECUTIVE STEPS 

2103 

1363 

PFIX, 

TAD 

PLEVW 

/GET CURRENT WORKING LEVEL 

2104 

4777 


JMS 

I CONVRP 

/AND CHANGE IT TO A PROBABILITY 

2105 

5774 


JMP 

I PNEXT 

/AND THAT'S IT 



/RE-ENTER PEST HERE TO GET NEW INPUT AND TYPE NEXT LEVEL 

2106 

4776 

PLOOP, 

JMS 

I PINPUT 

/GET A +1 OR -1 

2107 

3340 


DCA 

RESULT 

/SAVE THE RESULT 

2110 

1340 


TAD 

RESULT 

/NOW USE IT TO CALL PEST 

2111 

4750 


JMS 

I PESTP 


2112 

2163 


PLEVW / 

/WITH THIS PARAMETER LIST 

2113 

5316 


JMP 

PSTOP 

/IF RETURN IS MADE HERE, THAT'S IT 

2114 

7200 


CLA 

/ 

/IGNORE THE NEW LEVEL UNTIL LATER 

2115 

5232 


JMP 

PPRINT 

/PRINT CURRENT STATUS OF PEST 



/GO HERE IF 

STOPPING 

CRITERION MET 

2116 

4753 

PSTOP, 

JMS 

I RESETP 

/GO TO NEW LINE 

2117 

1351 

. 

TAD 

FINALP 


2120 

3340 


DCA 

RESULT 

/SETUP A POINTER TO A TEXT LIST 

2121 

1740 


TAD 

I RESULT 


2122 

7450 


SNA 



2123 

5327 


JMP 

FI NALL 


2124 

4752 


JMS 

I TYPEP 


2125 

2340 


ISZ 

RESULT 


2126 

5321 


JMP 

.-5 / 

/TYPE EACH NON-0 CHARACTER 

2127 

1363 

FI NALL, 

TAD 

PLEVW 


2130 

4746 


JMS 

I SPRINP 

/TYPE FINAL LEVEL 

2131 

4753 


JMS 

I RESETP 


2132 

7402 


HLT 

/ 

/THAT'S ALL (THIS OR THE NEXT 

2133 

7000 


NOP 

/ 

/LOCATION IS A GOOD PLACE FOR A 




/ 

/ 

/BREAKPOINT, OR TO RECYCLE FOR 




/ 

/ 

/RUN - TAKE YOUR CHOICE) 
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2134 

2233 

FIXALP, 

FIXALL 


2135 

7-671 

M71, 

-107 


2136 

0000 

LCOUNT, 

0 


2137 

0240 

SPACE, 

240 


2140 

0030 

RESULT, 

0 


2141 

0254 

C254, 

254 


2142 

7771 

M7, 

-7 


2143 

0056 

C56, 

56 


2144 

0304 

LETRD, 

304 


2145 

0255 

MINUS, 

255 


2146 

2310 

SPRINP, 

SPRINT 


2147 

1600 

WALDP, 

1600 


2150 

1637 

PESTP, 

1637 


2151 

2303 

FINALP, 

FINALT 


2152 

2215 

TYPEP, 

TYPE 


2153 

2223 

RESETP, 

RESET 


2154 

0260 

C260, 

260 


2155 

0000 

PGETW, 

0 / 

/ROUTINE TO USE WALD RESULT FOR PEST 

2156 

1340 


TAD RESULT 


2157 

5755 


JMP I PGETW 

/SIMPLE, ISN’T IT? 



/DUMMY ! 

SUBROUTINE TO 

CLEAR ACCUMULATOR 

2160 

0000 

CDUMMY, 

0 


2161 

7200 


CLA 


2162 

5760 


JMP I CDUMMY 




/WORKING PEST PARAMETER LIST SPACE 

2163 

0000 

PLEVW, 

0 / 

/CURRENT LEVEL 

2164 

0000 

PSTEPW, 

0 / 

/CURRENT STEP SIZE 

2165 

0000 

PMINW, 

0 / 

/MINIMUM STEP 

2166 

0000 

PMAXW, 

0 • / 

/MAXIMUM STEP 

2167 

0000 

PDIRW, 

0 / 

/CURRENT DIRECTION, -1 OR +1 FOR 




/ / 

/LAST STEP RESULTED FROM TOO MANY 




/ / 

/FAILURES OR SUCCESSES, RESPECTIVELY 

2170 

0000 

PDUBW, 

0 / 

/-I IF LAST STEP WAS A DOUBLING, 




/ / 

/0 OTHERWISE 

2171 

3000 

PCONW, 

0 / 

/TALLIES +1 FOR EACH SUBSEQUENT 


/ / /STEP AFTER THE FIRST IN ONE 

/ / /DIRECTION, STARTING FROM 0 AFTER A 

/ / /NON-DOUBLED STEP, FROM -1 IF THE 

/ / /NEW DIRECTION FOLLOWED A DOUBLED 

/ / /STEP 

PAUSE 



2173 

2174 

2175 

2176 

2177 


2200 

2201 

2202 

2203 


2204 

2205 

2206 
2207 


2210 

2211 

2212 

2213 

2214 

2215 

2216 
2217 
2220 
2221 
2222 

2223 

2224 

2225 

2226 
2227 

2230 

2231 

2232 


2233 

2234 

2235 


#2173 

/ALL THOSE LOCATIONS WHICH MIGHT NORMALLY BE CHANGED VIA 
/ODT ARE HERE PLACED TOGETHER, FROM 2173 TO 2207, FOR 
/THE CONVENIENCE OF THE PROGRAMMER: 


2106 

WNEXT, 

PLOOP 

/ 

2001 

PNEXT, 

WLOOP 

/ 

2245 

WINPUT, 

WACEPT 

/ 

2155 

PINPUT, 

PGETW 

/ 

2160 

CONVRP, 

CDUMMY 

/ 



/ 

/ 



/ 

/ 


/UNTIL CHANGED, THIS PROGRAM WILL 
/ALTERNATE WALD AND PEST TESTS, 
/TAKING WALD INPUT FROM THE TELETYPF 
/AND PASSING DECISIONS TO PEST 
/SPACE FOR ADDRESS OF ROUTINE TO 
/CONVERT CURRENT LEVEL TO A 
/PROBABILITY 


0144 
0020 
0032 
004 0 


DECIMAL 


/PEST STARTING PARAMETER LIST SPACE 


PEST0, 100 

16 

2 

32 


/ 

/ 

/ 

/ 


/STARTING LEVEL: EXAMPLE 100 
/STARTING STEP SIZE: EXAMPLE 16 
/MINIMUM STEP SIZE: EXAMPLE 2 
/MAXIMUM STEP SIZE: EXAMPLE 32 # 


/WALD PARAMETER LIST SPACE 



/REMEMBER: FRACTIONS 

DO NOT HAVE TO BE IN LOWEST TERMS 

0113 

WALD0 , 

75 / 

/TARGET NUMERATOR: EXAMPLE 75 

0144 


100 / 

/COMMON DENOMINATOR: EXAMPLE 100 

0226 

WDEVN, 

150 / 

/DEVIATION NUMERATOR: EXAMPLE 150 

0000 

WHIST, 

0 / 

/SPACE FOR WALD HISTORY 



OCTAL 


0000 

READ, 

0 


6031 


KSF 


5211 


JMP .-1 


6036 


KRB 


5610 


JMP I READ 


0000 

type. 

0 


6341 


TSF 


5216 


JMP .-1 


6046 


TLS 


7200 


CLA / 

/NOTE THAT THIS CLEARS ACCUMULATOR 

56 15 


JMP I TYPE 


0000 

RESET, 

0 


1231 


TAD CR 


4215 


JMS TYPE 


1232 


TAD LF 


4215 


JMS TYPE 


5623 


JMP I RESET 


0215 

CR, 

215 


0212 

LF, 

212 



/GRAND 

INITIALIZER 


0000 

FIXALL, 

0 


6046 


TLS / 

/START TELETYPE 

4644 


JMS I PINITP 

/INITIALIZE PEST 
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2236 

2200 


PEST0 / 

2237 

2163 


PLEVW / 

2240 

1206 


TAD WDEVN 

2241 

704 1 


CIA 

2242 

3237 


DCA WHIST 

2243 

5633 


JMP I FIXALL 

2244 

1715 

PINITP, 

1715 



/TYPEWRITER WALD INP 

2245 

0000 

WACEPT, 

0 

2246 

4210 


JMS READ 

2247 

7041 


CIA 

2250 

1260 


TAD A260 

2251 

7450 


SNA / 

2252 

5645 


JMP I WACEPT 

2253 

7031 


I AC / 

2254 

7640 


SZA CLA / 

2255 

5246 


JMP WACEPT+1 

2256 

7001 


I AC / 

2257 

5645 


JMP I WACEPT 

2260 

0263 

A263, 

260 

2261 

7777 

Ml, 

7777 


/PERMENANT PA 
/WORKING PARAM 


RAMETER LIST 
ETER LIST 


/INITIALIZE WALD 
/LEAVE 


ROUTINE 

/GET TYPEWRITER INPUT 

/WAS IT 260 (ASCII 0)? 
/IF SO, EXIT WITH 0 

/WAS IT 261 (ASCII 1)? 
/IF NOT, IGNORE IT 

/IF SO, RESTORE TO +1 
/AND EXIT 


2262 


2263 

2264 

2265 

2266 
2267 

2270 

2271 

2272 

2273 


2274 

2275 


2276 

2277 


0000 

4210 

1276 
7450 
5274 

1277 
7640 
5263 
7001 
5662 
7240 
5662 
7523 
0002 


/SUBROUTINE TO GET KEYBOARD INPUT FOR PEST. TYPE "+" FOR TOO 

K Ltf V f L L0WER> » R OR TOO LOW (MAKE LEVEL 

/HIGHER). THIS INPUT WILL NOT BE ECHOED. 

PACEPT, 0 


JMS READ /GET INPUT 

TAD MCMPAR /COMPARE TO 

SNA 


JMP TOOLOW 
TAD A2 
SZA CLA / 
JMP PACEPT+1 
TOOHI, I AC 

JMP I PACEPT 
TOOLOW, CLA CMA 

JMP I PACEPT 
MCMPAR, -255 
A2, 2 


/IF SO, TREAT ACCORDINGLY 

/WAS IT PLUS, THEN? 

/IF NOT, IGNORE IT 

/EXIT WITH +1 

/EXIT WITH -1 


2300 

0306 

FINALT, 306 

2301 

0311 

311 

2302- 

0316 

316 

2303 

0301 

301 

2304 

0314 

• 314 

2305 

0272 

272 

2306 

0240 

ASPACE, 240 

2307 

0000 

0 


/ASCII FOR "FINAL! " 


/"FINAL: " 


2310 

2311 

2312 


/DECIMAL PRINTING ROUTINE 
0000 SPRINT, 0 

3363 DCA NTEMP /SAVE NUMBER 

1362 TAD M3 
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2313 

3364 

DCA 

NHOLD 

2314 

1366 

TAD 

A 1000P 

2315 

3372 

DCA 

MPOINT 

2316 

3365 

DCA 

SPACER 

2317 

3373 

NLOOP1, DCA 

COUNTR 

2320 

7410 

SKP 

/ 

2321 

3363 

NL00P2, DCA 

NTEMP 

2322 

1772 

TAD 

I MPOINT 

2323 

7141 

CLL 

CIA / 

2324 

1363 

TAD 

NTEMP 

2325 

7420 

SNL 

/ 

2326 

5331 

JMP 

NOTFIT 

2327 

2373 

ISZ 

COUNTR 

2330 

5321 

JMP 

NL00P2 

2331 

7200 

NOTFIT, CLA 

/ 

2332 

2372 

ISZ 

MPOINT 

2333 

1373 

TAD 

COUNTR 

2334 

7450 

SNA 

/ 

2335 

5343 

JMP 

NZERO 

2336 

1260 

TAD 

A260 

2337 

4215 

JMS 

TYPE 

2340 

7001 

IAC 


2341 

3365 

DCA 

SPACER 

2342 

5350 

JMP 

NTALLY 

2343 

1365 

NZERO, TAD 

SPACER 

2344 

7650 

SNA 

CLA / 

2345 

1374 

TAD 

M20 / 

2346 

1260 

TAD 

A260 

2347 

4215 

JMS 

TYPE 

2350 

2364 

NTALLY, ISZ 

NHOLD 

2351 

5317 

JMP 

NLOOP1 

2352 

1363 

TAD 

NTEMP 

2353 

1260 

TAD 

A260 

2354 

4215 

JMS 

TYPE 

2355 

1306 

TAD 

ASPACE 

2356 

4215 

JMS 

TYPE 

2357 

1306 

TAD 

ASPACE 

2360 

4215 

JMS 

TYPE 

2361 

5710 

JMP 

I SPRINT 

2362 

7775 

M3, -3 


2363 

0000 

NTEMP, 0 


2364 

0000 

NHOLD, 0 


2365 

0000 

SPACER, 0 


2366 

2367 

A 1000P, A 1000 

2367 

1750 

A 1000, 1750 


2370 

0144 

144 


2371 

0012 

12 

/ 

2372 

0000 

MPOINT, 0 


2373 

0300 

COUNTR, 0 


2374 

7760 

M20, -20 



PAUSE 


/SAVE A -3 COUNTER 

/SETUP POINTER TO LIST OF MINUENDS 

/PRINT 3LANKS FOR LEADING 0S 

/SETUP 0 COUNT 

/DON’T DESTROY THE NUMBER 

/SAVE THE NUMBER 

/GET MINUEND 

/CLEAR LINK AND NEGATE 

/ADD CURRENT VALUE OF NUMBER 

/OVERFLOW? 

/IF NOT, MINUEND WON’T FIT 

/IF SO, TALLY 1 

/AND TRY AGAIN 

/LEAVE NTEMP AS IS 

/SETUP NEXT MINUEND 

/NOW TRY COUNTR 

/IF NONZERO, ALWAYS PRINT 

/IF ZERO, TEST FOR SPACE OR ”0" 

/ADD DIGIT VALUE TO 260 

/OUTPUT IT 

/MAKE SPACER NON-0 
/TEST FOR 3 *D DIGIT 
/GET SPACER 
/AND, IF IT’S 0, 

/ADD -20 TO 260, MAKING 240 

/OUTPUT IT 
/3 DONE YET? 

/IF NOT, DO MORE 

/IF SO, PRINT LAST DIGIT 

/TYPE IT 


/TYPE 2 SPACES J 

/EXIT 


/1000, 100, 10 IN DECIMAL 
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*2400 

/ROUTINE TO GENERATE A +1 WITH RATIONAL PROBABILITY 
/EQUAL TO PNUtl/PDEN, TREATED AS 12-31T NON-NEGATIVE INTEGERS 
/THIS ROUTINE CAN 3E CALLED BY PLACING ITS ADDRESS INTO 
/LOCATION WIMPUT. THE VALUES PNUM AND PDEN ARE TO BE FILLED 
/BY A ROUTINE WHOSE ADDRESS IS PLACED INTO CONVRP, AND 
/WHICH IS CALLED AFTER EACH CHANGE OF PEST LEVEL. 


2400 

0000 

PNUM, 

0 

/ 

/NUMERATOR AND DENOMINATOR OF 

2401 

0003 

PDEN, 

0 

/ 

/PROBABILITY OF GENERATING A +1 

2402 

0000 

PR0B1, 

0 



2433 

7200 


CLA 

/ 

/TO BE SAFE 

2404 

1201 


TAD 

PDEN 

/GET DENOMINATOR 

2405 

7450 


SNA 



2406 

7402 


HLT 

/ 

/MAKE SURE IT IS NON-0 

2437 

4215 


JMS 

RANDOM 

/CHOOSE FROM 1 TO IT AT RANDOM 

2410 

7141 


CLL 

CIA / 

/NEGATE AND CLEAR LINK 

2411 

1200 


TAD 

PNUM 

/ADD NUMERATOR - IF LINK CHANGES, 




/ 

/ 

/THEN NEGATED RANDOM RESULT WAS 




/ 

/ 

/NO MORE NEGATIVE THAN NUMERATOR 




/ 

/ 

/IS POSITIVE, SO WE GENERATE A +1 




/ 

/ 

/FOR A SUCCESS 

2412 

7630 


SZL 

CLA 


2413 

7001 


IAC 



2414 

5602 


JMP 

I PR0B1 

/EXIT WITH +1 OR WITH 3 



/ROUTINE TO 

CHOOSE A PSEUDO-RANDOM NUMBER FROM THE RANGE 



/I TO N 

, WHEN ROUTINE 

IS ENTERED WITH N IN THE ACCUMULATOR 



/ROUTINE EXITS WITH RANDOM NUM3ER I N ACCUMULATOR 




/IF 

N>1, WORK 

WITH RANGE 0 TO (N-l) FOR NOW 

2415 

0000 

RANDOM, 

0 



24 IS 

7450 


SNA 



2417 

5615 


JMP 

I RANDOM 


2420 

7041 


CIA 



2421 

7040 


CMA 



2422 

7450 


SNA • 



2423 

5254 


JMP 

INCSEL 





/FIND LEFTMOST 

1 BIT AND CREATE MASK OF 1 BITS FROM 




/THERE TO THE 

RIGHT END OF THE ACCUMULATOR 

2424 

3256 


DC A 

LIMIT 


2425 

3257 


DCA 

BITTA3 


2426 

1256 


TAD 

LIMIT 


2427 

7110 

SHIFT, 

CLL 

RAR 


2430 

2257 


ISZ 

BITTAB 


2431 

7440 


SZA 



2432 

5227 


JMP 

SHIFT 


2433 

1257 


TAD 

3ITTAB 


2434 

7041 


CIA 



2435 

3257 


DCA 

BITTAB 


2436 

7120 

MASKER, 

STL 



2437 

7004 


RAL 





2440 

2257 


ISZ BITTA8 


2441 

5236 


JMP MASKER 


2442 

3260 


DCA MASC 





/GET 12-BIT RANDOM NUMBER, MASK 

TO APPROPRIATE 




/NUMBER OF BITS, KEEP IT IF IT 

IS NOT TOO LARGE, 




/OTHERWISE TRY AGAIN 


2443 

4261 

RNGET, 

JMS FLIP 


2444 

0260 


AND MASC 


2445 

3257 


DCA CHOICE 


2446 

1256 


TAD LIMIT 


2447 

7140 


CMA CLL 


2450 

1257 


TAD CHOICE 


2451 

7630 


SZL CLA 


2452 

5243 


JMP RNGET 


2453 

1257 


TAD CHOICE 





/ADD +1 TO RESTORE TO SELECTION 

FROM RANGE 1 TO N 

2454 

7001 

INCSEL, 

IA C 

2455 

5615 


JMP I RANDOM 


2456 

0000 

LIMIT, 

0 




BITTABr 

• 


2457 

0000 

CHOICE, 

0 


2460 

0000 

MASC, 

0 



/SUBROUTINE TO GENERATE 12-BIT PSEUDO-RANDOM NUMBER VIA THE 
/RECURSION X(N+1)=X(N)+X(N-21), MODULUS (2tl2) 


2461 

0000 

FLIP, 

0 

2462 

7200 


CLA 

2463 

1302 


TAD FLPT1 

2464 

3301 


DCA FLPT 

2465 

1302 


TAD FLPT1 

2466 

1303 


TAD FLBOT 

2467 

7640 


SZA CLA 

2470 

5273 


JMP FLOK 

2471 

1304 


TAD FLTOP 

2472 

3302 


DCA FLPT 1 

2473 

2302 

FLOK, 

ISZ FLPT1 

2474 

1701 


TAD I FLPT 

2475 

1702 


TAD I FLPT 1 

2476 

3702 


DCA I FLPT 1 

2477 

1702 


TAD I FLPT 1 

2500 

5661 


JMP I FLIP 

2501 

2505 

FLPT, 

LTOP 

2502 

2506 

FLPT1, 

LTOP+1 

2503 

5246 

FLBOT, 

-LBOT 

2504 

2504 

FLTOP, 

LTOP-1 

2505 

2315 

LTOP, 

2315 

2506 

2303 


2303 

2537 

7553 


7553 

2510 

3676 


3676 

2511 

0655 


0655 

2512 

4331 


4331 

2513 

5255 


5255 



2514 

3135 

3135 

2515 

3157 

3157 

2516 

5704 

5704 

2517 

3072 

3072 

2520 

5223 

5223 

2521 

5667 

5667 

2522 

7262 

7262 

2523 

3566 

3566 

2524 

6102 

6102 

2525 

4547 

4547 

2526 

3535 

3505 

2527 

6475 

6475 

2530 

1047 

1047 

2531 

755 1 

7551 

2532 

1132 LBOT, 

1 132 
PAUSE 
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4 


/THIS IS AN EXAMPLE OF A SUBROUTINE TO CONVERT THE CURRENT 
/LEVEL TO A PROBABILITY, FOR SIMULATING THE ACTION OF PEST. 
/IT CREATES A VERY CRUDE APPROXIMATION TO A CUMULATIVE 
/NORMAL OGIVE BY ADJOINING PARTS OF TV/O PARABOLAS, ONE 
/IN NORMAL ORIENTATION, ONE INVERTED. THIS APPROXIMATION 
/GENERATES PROBABILITIES FROM A LOCATION PARAMETER, C, WHICH 
/DETERMINES THE LEVEL AT WHICH PR03(SUCCESS) = .75. FOR 
/LEVELS BETWEEN C-S4 AND C, THE CURVE IS A SEGMENT OF THE 
/RIGHT HALF OF AN UPWARD-OPENING PARABOLA WHICH IS TANGENT 
/TO THE LINE PROS(SUCCESS>=.5 AT THE POINT (C-54,.5), 

/AND WHICH MEETS THE LEFT HALF OF A CONGRUENT BUT INVERTED 
/PARABOLA SEGMENT AT (C,.75). THIS INVERTED SEGMENT IS 
/TANGENT TO THE LINE PROB(SUCCESS>=1.0 AT THE POINT 
/(C+64,1.0). BELOW OR ABOVE THE RANGE fC-G4,C+64] THE 
/PROBABILITY OF A SUCCESS IS .5 OR 1.0, RESPECTIVELY. 

/THUS THE GRAPH IS SIMPLY .5 + THE INTEGRAL OF AN ISOCELES 
/TRIANGLE WITH BASE 128, AREA .5, CENTERED AT C. 

/HERE IS A GRAPH OF THE FUNCTION, AS PRODUCED BY FOCAL: 


/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

/ 


1.00 

0.95 

0.90 

0.85 

0.80 

0.75 

0.70 

0.55 

0.50 

0.55 

0.50 

0.45 

P 


• • • •+(C+54,1.) # 

• • • 

• • • 


+(C,.75) 


(C-64,.5) 


C-54 C 

LEVEL 


* 



C+54 


/WITHIN THE RANGE [-54,541, THE PROBABILITY FUNCTION IS 
/DEFINED AS FOLLOWS: 

/IF X<=C, P(SUCCESS) = .5 + [(X-(C-64))T2]/(2f14), AND 
/IF X>=C, P(SUCCESS) = 1. - [(X-(C+64))t21/(2tl4). 
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/FOR COMPUTATIONAL PURPOSES EACH OF THESE CAN BE EXPRESSED 
/AS A FRACTION WITH DENOMINATOR 2tll (4000 OCTAL), AND 
/THE FORMULAS FOR THE NUMERATORS REDUCE TO 

/IF X<=C, NUMERATOR = 2tl0 + r (X-(C-64))t2]/8, AMD 
/IF X>=C, NUMERATOR = 2t l l - [(X-(C-64))t2]/?♦ 

*2600 




2600 

0000 

CURVE, 

0 

/ 

2601 

3262 


DCA 

CLEVEL 

2602 

1266 


TAD 

C 

2603 

1263 


TAD 

A 64 

2604 

7141 


CLL 

CIA / 

2605 

1262 


TAD 

CLEVEL 

2606 

7470 


SZL 

SNA / 

2607 

5224 


JMP 

SET 10 

2610 

1263 


TAD 

A64 / 

261 1 

7430 


SZL 

/ 

2612 

5226 


JMP 

UPHALF 

2613 

1263 


TAD 

A 64 / 

2614 

7430 


SZL 

/ 

2615 

5234 


JMP 

LOHALF 

2616 

7332 

SET05, 

CLA 

STL RTR 

2617 

3664 

EXIT, 

DCA 

I NUMPTR 

2620 

1256 


TAD 

A 4000 

2621 

3665 


DCA 

I DENPTR 

2622 

7000 


NOP 

/ 




/ 

/ 

2623 

5600 


JMP 

I CURVE 

2624 

7330 

SET10, 

CLA 

STL RAR 

2625 

5217 


JMP 

EXIT 

2626 

704 1 

UPHALF, 

CIA 

/ 




/ 

/ 

2627 

1263 


TAD 

A 64 / 




/ 

/ 




/ 

/ 

2630 

4237 


JMS 

SQUARE 

2631 

7041 

. 

CIA 

/ 

2632 

1256 


TAD 

A 4000 

2633 

5217 


JMP 

EXIT 

2634 

4237 

LOHALF, 

JMS 

SQUARE 




/ 

/ 

2635 

1257 


TAD 

A2000 

2636 

5217 


JMP 

EXIT 

2637 

0000 

SQUARE,' 

0 

/ 




/ 

/ 




/ 

/ 




/ 

/ 

2640 

3262 


DCA 

CLEVEL 

2641 

1262 


TAD 

CLEVEL 

2642 

7141 


CLL 

CIA / 

2643 

7450 


SNA 

/ 


/ENTER WITH LEVEL IN ACCUMULATOR 
/SAVE LEVEL 


/CLEAR LINK TO RECEIVE OVERFLOW 
/LEVEL >= C+64 ? 

/YES, IF LINK CHANGES OR ACC = 0 
/SET PR03 TO 1.0 IF SO 
/C < LEVEL < C+64 ? 

/YES, IF LINK CHANGES NOW 
/IF SO, IN UPPER HALF OF CURVE 
/C-64 < LEVEL <= C? 

/YES, IF LINK CHANGES NOW 
/IF SO, IN LOWER HALF OF CURVE 
/CREATE UNCONDITIONAL 2300 NUMERATOR 
/STORE NUMERATOR 

/USE 4000 DENOMINATOR ALWAYS 
/REPLACE WITH JMS WRITIT TO PRINT 
/THE NEW LEVEL IN DECIMAL 
/RETURN 

/CLEAR ACCUMULATOR AND CREATE 4000 

/FOR PROBABILITY OF ,5 

/ACCUMULATOR CONTAINED LEVEL - C, SO 

/CHANGE TO C - LEVEL, THEN 

/ADD 64 FOR (C+64) - LEVEL, WHICH IS 

/POSITIVE BUT HAS THE SAME SQUARE AS 

/LEVEL - (C+64) 

/GET ITS SQUARE, DIVIDED BY 8 
/NEGATE IT 
/ADD 2t11 
/AND GO 

/ACCUMULATOR HAD L-(C-64), SO JUST 
/GET SQUARE DIVIDED 3Y 8 
/ADD TO 2t10 
/AND GO 

/ROUTINE TO GET 13-BIT SQUARE 
/OF NUMBER <= 64, AND TO DIVIDE 
/THE RESULT BY 8, WITH APPROPRIATE 
/ROUNDING 
/SAVE NUMBER 

/TRY TO CLEAR LINK, AND NEGATE 
/BUT IF NUMBER WAS 0, 
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2644 

5637 


JMP I SQUARE 

/LEAVE WITH A 0 

2645 

3261 


DCA COUNTX 

/OTHERWISE, SAVE A COUNTER 

2646 

1262 


TAD CLEVEL 

/NOW MULTIPLY BY ITSELF 

2647 

2261 


ISZ COUNTX 


2650 

5246 


JMP .-2 


2651 

1260 


TAD A4 / 

/ADD 4 WHEN DONE 

2652 

7010 


RAR / 

'/SAVE LINK THE FIRST TIME 

2653 

7110 


CLL RAR 


2654 

7110 


CLL RAR / 

/THEN DIVIDE BY 4 MORE 

2655 

5637 


JMP I SQUARE 

/YIELDING ROUNDED 8-BIT SQUARE 

2656 

4000 

A 4000, 

4000 


2657 

2000 

A2000, 

2000 


2660 

0004 

A 4, 

4 


2661 

0000 

COUNTX, 

0 


2662 

0000 

CLEVEL, 

0 


2663 

0100 

A 64 , 

100 


2664 

2400 

NUMPTR, 

PNUM 


2665 

2401 

DENPTR, 

PDEN 


2666 

0400 

c, 

400 / 

/C MAY 3E ANY INTEGER IN THE RANGL 




/ / 

PAUSE 

/101 TO 7677, OCTAL 
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/OPTIONAL SUBROUTINE TO PRINT THE CURRENT PROPORTION AS A 
/FRACTION WITH DENOMINATOR 1000 

/THIS ROUTINE BEGINS BY CALCULATING T, THE NUMBER OF 
/THOUSANDTHS EQUAL TO THE NEW PROPORTION. FROM THE FORMULA 

/ T = INTEGER PART OF ((NUMERATOR * 1000 / 2048) + .5) 

/WE CAN DERIVE THE COMPUTATIONAL FORMULA 

/ T = [ ((2048*(NUMERATOR+1> - 48*NUMERAT0R) / 4096)]. 

/THIS LATTER FORMULA, OF COURSE, IS VERY EASY TO IMPLEMENT 


/ON A 12-BIT COMPUTER. 
/DECIMAL, FOLLOWED BY 



2667 

0000 

WRITIT, 

0 



2670 

1664 


TAD 

I NUMPTR 

• 

2671 

7001 


IAC 

/ 


2672 

7110 


CLL 

RAR 


2673 

3332 


DCA 

TOP 


2674 

7010 


RAR 



2675 

3333 


DCA 

BOT / 


2676 

1334 


TAD 

M48 


2677 

3261 


DCA 

COUNTX 


2700 

1664 


TAD 

I NUMPTR 


2701 

7041 


CIA 



2702 

3331 


DCA 

COPNUM 


2703 

1333 

LOOPW, 

TAD 

BOT / 


2704 

1331 


TAD 

COPNUM 


2705 

3333 


DCA 

BOT 


2706 

7420 


SNL 



2707 

7040 


CMA 

/ 


2710 

1332 


TAD 

TOP 


2711 

3332 


DCA 

TOP / 


2712 

7100 


CLL 



2713 

2261 


ISZ 

COUNTX 

• 

2714 

5303 


JMP 

LOOPW 


2715 

1335 


TAD 

OPAREN 


2716 

4747 

' 

JMS 

I TYPE2 


2717 

1332 


TAD 

TOP 


2720 

4750 


JMS 

I SPRIN2 


2721 

1336 


TAD 

CLOSEP 


2722 

3261 


DCA 

COUNTX 


2723 

1661 


TAD 

I COUNTX 


2724 

7450 


SNA 



2725 

5667 


JMP 

I WRITIT 


2726 

4747 


JMS 

I TYPE2 


2727 

2261 


ISZ 

COUNTX 


2730 

5323 


JMP 

.-5 


2731 

0000 

COPNUM, 

0 



2732 

0000 

TOP, 

0 



2733 

0000 

30T, 

0 



THE VALUE OF T IS THEN PRINTED IN 

"/ 1000 ". 

/ADD +1 TO NUMERATOR AND 

/MULTIPLY BY 2048, DOUBLE PRECISION 
/PREPARE TO SUBTRACT 48 TIMES 

/CREATE NEGATIVE OF NUMERATOR 
/ADD NEGATIVE TO BOTTOM 48 TIMES, 

/SUBTRACTING 1 FROM THE TOP 
/EACH TIME THE LINK DOESN’T CARRY 

/PRINT "(" 

/PRINT TOP 

/PRINT "/ 1000)" UNTIL DONE 
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2734 

7720 

M48, 

-60 

2735 

0250 

OPAREN, 

250 

2735 

2737 

CLOSEP, 

• +1 

2737 

0257 


257 

2743 

0240 


240 

2741 

0261 


261 

2742 

0260 


260 

2743 

0260 


260 

2744 

0260 


260 

2745 

0251 


251 

2746 

0030 


0 

2747 

2215 

TYPE2, 

TYPE 

2750 

2310 

SPRIN2, 

SPRINT 



ASPACE 

2306 

A 1 033 

2367 

A1003P 

2366 

A2 

2277 

A2000 

2657 

A260 

2263 

A 4 

2663 

A4000 

2656 

A64 

2663 

BITTAB 

2457 

BOT 

2733 

C 

2666 

CDUMMY 

2160 

CHOICE 

2457 

CLEVEL 

2662 

CLOSEP 

2736 

CONVRP 

2177 

COPNIJM 

2731 

COUNTR 

23 73 

COUNTX 

2661' 

CR 

2231 

CURVE 

2600 

C254 

2141 

C260 

2154 

C56 

2143 

DENPTR 

2665 

EXIT 

2617 

FI NALL 

2127 

FINALP 

2151 

FI NALT 

2303 

FI XALL 

2233 

FIX ALP 

2134 

' FLBOT 

2503 

FLIP 

2461 

FLOK 

2473 

FLPT 

2501 

FLPT1 

2502 

FLTOP 

2504 

INCSEL 

2454 

LBOT 

2532 

LCOUNT 

2136 

LETRD 

2144 

LF ’ 

2232 

LIMIT 

2456 

LOHALF 

2634 

LOOPW 

2703 

LTOP 

2505 

MASC 

2460 

MASKER 

2436 

MCMPAR 

2276 

MINUS 

2145 

MPOINT 

2372 

Ml 

2261 



M2 0 

2374 

>13 

2362 

M48 

2734 

M7 

2142 

M71 

2135 

NEGCT 

2071 

NHOLD 

2364 

NLOOP1 

2317 

HL00P2 

2321 

NOTFIT 

2331 

NTALLY 

2350 

NTEMP 

2363 

NUMPTR 

2564 

nzero 

2343 

OPAREN 

2735 

PACEPT 

2262 

PCONW 

2171 

PDEN 

2401 

PDIRW 

2167 

PDUBW 

2170 

PESTGO 

2030 

PESTP 

2150 

PEST0 

2200 

PFIX 

2103 

PGETW 

2155 

PINITP 

2244 

PINPUT 

2176 

PLEVW 

2163 

PLOOP 

2106 

PMAXW 

2166 

Pill NW 

2165 

PNEXT 

2174 

PNUM 

2400 

POSCT 

2077 

PPRINT 

2032 

PRINTC 

2045 

PRINTL 

2051 

PR OBI 

2402 

PSTEPW 

2164 

PSTOP 

2116 

RANDOM 

2415 

READ 

2210 

RESET 

2223 

RESETP 

2153 

RESULT 

2140 

RNGET 

2443 

SET05 

2616 

SET 10 

2624 

SHIFT 

2427 

SPACE 

2137 

SPACER 

2355 

SPRINP 

2146 

SPRINT 

2310 

SPRIN2 

2750 



SQUARE 

2637 

TOOHI 

2272 

TOOLOW 

2274 

TOP 

2732 

TYPE 

2215 

TYPEP 

2152 

TYPE2 

2747 

UPHALF 

2626 

WACEPT 

2245 

WAGAIN 

2005 

WALDGO 

2003 

WALDP 

2147 

WALDO 

2204 

WDEVN 

2206 

WHIST 

2207 

WINPUT 

2175 

WLINE 

2002 

WLOOP 

2001 

WNEXT 

2173 

WRITIT 

2667 
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PEST DEMONSTRATOR -- ANNOTATED DEMONSTRATION 


THE TYPED MATERIAL OCCUPYING, ROUGHLY, THE LEFT HALF OF THESE 
PAGES IS THE OUTPUT OF ODT AND THE DEMONSTRATOR. ECHOED OPERATOR INPUT 
IS UNDERLINED. ALL PARENTHETICAL COMMENTS ARE ANNOTATIONS ADDED LATER, 
AS ARE THESE VARIOUS INTRODUCTORY AND INTERSPERSED PARAGRAPHS. 


TO START THESE DEMONSTRATIONS, THE PEST PACKAGE (AT 1S00), THE PEST 
DEMONSTRATOR PACKAGE, AND ODT LOW WERE LOADED USING THE BIN LOADER. FOR 
THE FIRST DEMONSTRATION, THE PROGRAM WAS USED AS LOADED, WITH WALD 
GETTING STRINGS OF "1”S AND "0"S FROM THE OPERATOR, MAKING A DECISION 
ABOUT EACH STRING, AND PASSING THE DECISION TO THE DEMONSTRATOR TO BE 
PRINTED AND TO PEST FOR A NEW LEVEL TO BE CHOSEN AND PRINTED. 

THE PROGRAM WAS STARTED AT 2030 BY ODT. 


2030G 


(START VIA ODT) 


100 ND 0 (THE DEMOSTRATOR TYPES THE CURRENT LEVEL, 

WHETHER OR NOT THE STEP WAS DOUBLED, AND THE 
VALUE OF PEST’S INTERNAL COUNTER TO DECIDE 
WHEN TO BEGIN DOUBLING THE STEP SIZE.) 


000 * 


+ 16 

116 

ND 

1 

000— 

+ 32 

148 

D 

2 

000- 
+ 32 

180 

D 

3 

1111111+ 



- 16 

164 

ND - 

1 

1111111+ 



--— 

148 

ND 

0 

100- 
+ 8 

156 

ND 

0 

1111111+ 

4 

152 

ND 

0 

1111111+ 

4 

148 

ND 

1 

000- 
+ 2 

150 

ND 

0 

1111111+ 

FINAL: 

150 




(THE OPERATOR TYPES "000" FOR THREE CONSECUTIVE 
FAILURES, AND WALD RESPONDS WITH "TOO LOW".) 

(SO PEST RAISES THE LEVEL BY 16, TO 116) 


(NOTE THE DOUBLED STEP.) 


(NO FURTHER DOUBLING, AS THE MAXIMUM IS 32.) 

(OPERATOR TYPES A STRING OF SUCCESSES, TO CAUSE 
A REVERSAL OF DIRECTION. STEP SIZE DROPS TO 16.) 


(OPERATOR CAUSES CONVERGENCE AT 150, BECAUSE 
THE MINIMUM STEP ALLOWED WAS SET TO 2.) 
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NOW ODT IS USED TO MAKE WALD SIMPLY RETURN TO ITSELF WHEN DONE, SO 
THAT THE OPERATOR CAN SEE WHAT DECISION WALD MAKES ON SOME OF THE 
POSSIBLE INPUT STRINGS. THE OPERATOR TYPES A STRING OF IS AND 0S, 

AND WALD RESPONDS "+" OR WALD IS STILL SET FOR TARGET LEVEL OF 757., 

DEVIATION LIMIT OF 1.5 TRIALS FROM EXPECTATION. REMEMBER THAT WALD 
IGNORES THE FIRST TRIAL OF ANY STRING. 


2173/ 2136 2001 (ODT USED TO CHANGE WALD EXIT POINTER.) 

2000G 


000 - 


01000 “ 
011000 - 


0111111 + 


00111111111 + 


(2 FAILURES IN 2 COUNTED TRIALS, WHILE ONLY 
.5 FAILURES WERE EXPECTED. DIFFERENCE IS 1.5. SO 
DECIDE THAT SUCCESS RATE IS TOO LOW.) 

(3 FAILURES IN 4 TRIALS, EXPECTING 1.0, TOO LOW) 

(3 FAILURES, 5 TRIALS, EXPECTING 1.25, TOO LOW) 

(6 SUCCESS, 6 TRIALS, EXPECTING 4.5, 

DIFFERENCE EXACTLY 1.5, TOO HIGH) 

(AFTER FIRST FAILURE, NEED 9 SUCCESS TO GET 
9 SUCCESSES OUT OF 10, EXPECTING 7.5, TOO HIGH) 


NOW ODT IS USED TO CHANGE THE DEVIATION LIMIT BACK TO 1.0 
2206/ 0226 144 . (144 WAS COMMON DENOMINATOR ALREADY) 

2000G 


000 - 

0100 - 


01100 - 


01111 + 

011101110111011101110111011101111 + 

(NOTE THAT 3 SUCCESSES PER FAILURE IS RIGHT 
ON TARGET, AND A DECISION COULD BE POSTPONED 
INDEFINITELY.) 
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NOW ODT IS USED TO RESET THE WALD EXIT POINTER AND TO CAUSE PEST TO 
LOOP 3ACK TO ITSELF, SO THAT THE OPERATOR CAN ENTER M +" OR TO 
PEST TO SEE THE EFFECTS OF VARIOUS SEQUENCES OF SUCH ON THE 
CHANGES OF PEST LEVEL. THE "+" OR IS NON-ECHOING, BUT THF SIGN 
ON THE CHANGE OF LEVEL IS ALWAYS THE OPPOSITE OF THE INPUT SIGN, 

SINCE "+" INPUT MEANS LEVEL TOO HIGH, SO THE NEXT LEVEL IS LOWER, ETC. 


2173/2001 2106 
2174/ 2001 2106 

2330G 


(RESTORE WALD 
(PROGRAM WILL 
TESTS) 


EXIT POINTER) 
SIMPLY LOOP FOR 


CONSECUTIVE PEST 




100 

ND 

0 

+ 

16 

116 

ND 

1 

+ 

32 

148 

D 

2 

- 

16 

132 

ND - 

1 

- 

16 

115 

ND 

0 

- 

16 

100 

ND 

1 

- 

32 

68 

D 

2 


- 

32 

36 

D 

3 

mm 

32 

4 

D 

4 

+ 

16 

20 

ND - 

1 

- 

8 

12 

ND 

0 

- 

8 

4 

ND 

1 

- 

16 

4084 

D 

2 


(OPERATOR HAD TYPED 


FOR LEVEL TOO LOW) 


(NOTICE THAT IT IS NOT UNTIL THE 4’TH 
CONSECUTIVE STEP AFTER A DOUBLING THAT THE 
STEP SIZE IS DOUBLED.) 


(TWO THINGS TO NOTICE HERE: FIRST, SINCE THE 
LAST REVERSAL DID NOT FOLLOW A DOUBLED STEP, THE 
THIRD CONSECUTIVE STEP IN THIS DIRECTION WAS 
DOUBLED. SECOND, THE LEVEL UNDERFLOWED, AN 
EVENT WHICH THE DEMONSTRATOR IGNORES.) 


+ 8 4092 ND 

+ 8 4 ND 

+ 8 12 ND 

-4 8 ND 

+ 2 10 ND 

FINAL: 10 


1 

0 

1 

0 

0 


(CONVERGENCE 


TO LEVEL OF 10) 
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NOW THE PROGRAM IS SET TO MAKE USE OF THE AUTOMATIC INPUT SIMULATOR 
WHICH GENERATES SUCCESSES AND FAILURES ACCORDING TO A SIMPLE * 

PSYCHOMETRIC FUNCTION, IMPLEMENTED 3Y SUBROUTINE "CURVE". THF DEVIATION 
LIMIT REMAINS AT 1.0 TRIALS. ' ViHllUN 


2174/ 2106 2001 

2175 /2245 2402 

2176 /2262 2155 

2177 /2313 


2030G 


(AFTER PEST, GO 3ACK TO WALD) 

(WALD GETS ITS INPUT FROM ROUTINE "PR031") 

(PEST GETS ITS INPUT FROM WALD) 

(THE ROUTINE "CURVE" IS TO 3E CALLED 
AFTER EACH CHANGE OF LEVEL, TO COMPUTE THE 
NEW PROBABILITY OF SUCCESS) 

(FROM THIS POINT ON, THE OPERATOR ENTERS 
NOTHING. THE ROUTINES "CURVE” AND "PR031" SUPPLY 
WALD WITH ITS SUCCESSES AND FAILURES. THE 
CALIBRATION OF CURVE IS PRESENTLY SET SO THAT A 
LEVEL OF 256 GENERATES THE TARGET 752 CORRECT. 

A LEVEL OF 192 OR LOWER GENERATES SUCCESSES 
WITH PROBABILITY 5 07., 323 OR HIGHER GENERATES 
1W7. SUCCESSES.) 



100 

ND 

0 

0130- 
+ 16 

116 

ND 

1 

100- 
+ 32 

148 

D 

2 

1313- 
+ 32 

180 

D 

3 

00111011132- 
+ 32 212 

D 

4 

11111-4- 
- 16 

196 

ND - 

1 

10111101100- 
+ 8 204 

ND 

0 

000- 
+ 8 

212 

ND 

1 

0011 3- 
+ 16 

228 

D 

2 

11111+ 

8 

220 

ND - 

1 

1103- 
+ 4 

224 

ND 

0 


0 (CONTINUED ON NEXT PAGE) 
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1313- 
+ 4 

228 

110111010- 

+ 8 

236 

11101131131 

+ 16 

252 

130- 
+ 32 

284 

11111+ 

- 16 

268 


D 

D 


1 

2 

3 

4 


0 


1110111101111 + 

- IS 252 ND 

101111313- 

+ 8 260 ND 


101 111101110110111011111 1 + 


4 

256 

ND 

0 

111101111+ 

4 252 

ND 

1 

000* 

+ 2 

254 

ND 

0 

1100- 

+ 2 

256 

ND 

1 


1101111011101111 1 + 

FINAL: 256 (SHEER LUCK) 

NOW THREE CHANGES ARE MADE: 

1. THE DEVIATION LIMIT IS RAISED TO 1.5 AGAIN. 

2. THE TARGET PERCENTAGE IS CHANGED TO 85%. 

3. THE "CURVE" ROUTINE PRINTS THE SUCCESS PROBABILITY CORRESPONDING TO 
EACH NEW LEVEL THAT IT SETS. THIS IS PRINTED IN PARENTHESES — THE 
FRACTIONS WITH 1000 IN THE DENOMINATOR ARE NOT ANNOTATIONS ADDED 
LATER. 

A 

2234/ 0113 125 
2206/ 0144 226 
2622/ 7000 4267 

2033G 

100 ND 


(CHANGE TARGET NUMERATOR FROM 75 TO 85) 
(CHANGE DEVIATION LIMIT TO 1.5) 

(IMPLEMENT PROBABILITY-PRINTING ROUTINE) 

(AGAIN, OPERATOR IS "HANDS-OFF" AFTER THIS) 

0 ( 500 / 1000) (CONTINUED ON NEXT PAGE) 
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001111010 - 


+ 16 116 

ND 

1 

( 500 

/ 

1000) 

1310- 

+32 148 

D 

2 

( 500 

/ 

1000) 

110103- 
+ 32 180 

D 

3 

( 500 

/ 

1000) 

001100- 
+ 32 212 

D 

4 

( 524 

/ 

1000) 

0131010- 
+ 32 244 

D 

5 

( 665 

/ 

1000) 

1011010- 
+ 32 276 

D 

6 

( 882 

/ 

1000) 

11111111111+ 
- IS 260 

ND - 

1 

( 780 

/ 

1000) 

0010- 

+ 8 268 

ND 

0 

( 835 

/ 

1000) 

11 1 101111111011111110111101 1 1 1 1 1 1 1 1 10! 101110010- 
+ 8 276 ND 1 ( 882 / 1000) 

11311111100- 
+16 292 

D 

2 

( 952 

/ 

1000) 

11011111111111111011 11 11 1 + 

8 284 ND - 1 ( 921 

/ 

1000) 

1011111111111110111101111111111+ 

8 276 ND 0 ( 882 / 

1000) 

1111111110111101111111111+ 

8 268 .ND 1 ( 835 

/ 

1300) 

11011111011111110011110- 
+ 4 272 ND 0 

( 859 

/ 

1300) 

01111111111+ 
2 270 

ND 

0 

< 847 

/ 

1000) 


11011111110110111111111111111011111111 + 

2 268 ND 1 C 835 / 1000) 

011011 111 11 11 11 101 111001 111 101 1111 11 11 111 11 1011011-30 10- 
FINAL: 268 


NOT TOO BAD --THE SUCCESS PROPORTION AT CONVERGENCE IS 83.5%, COMPARED 
TO A TARGET OF 85%. TO GET MORE INFORMATION ABOUT THE DISTRIBUTION 
OF FINAL LEVELS AND THEIR CORRESPONDING SUCCESS PROPORTIONS, THE MOST 
EFFICIENT USE OF THIS SIMULATOR WOULD BE TO DELETE MOST OF THE 
TYPED OUTPUT, RETAINING ONLY THE TERMINAL LEVELS AND, POSSIBLY, 
INFORMATION ABOUT THE NUMBER OF TRIALS NEEDED TO ACHIEVE CONVERGENCE. 
BUT HERE ENDS OUR PART OF THE DEMONSTRATION. 
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